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We calculated the first two moments of the light-cone distribution amplitudes for the pseudoscalar 
mesons (n and K) and the longitudinally polarised vector mesons (p, K* and <p) as part of 
the UKQCD and RBC collaborations' Nf = 2 + 1 domain- wall fermion phenomenology pro- 
gramme. These quantities were obtained with a good precision and, in particular, the expected 
effects of SU(3) -flavour symmetry breaking were observed. Operators were renormalised non- 
perturbatively and extrapolations to the physical point were made, guided by leading order chiral 
perturbation theory. The main results presented are for two volumes, 16 3 x 32 and 24 3 x 64, with a 
common lattice spacing. Preliminary results for a lattice with a finer lattice spacing, 32 3 x 64, are 
discussed and a first look is taken at the use of twisted boundary conditions to extract distribution 
amplitudes. 
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1. Introduction 

These proceedings will give our current values for the lowest moments of the leading twist 
meson distribution amplitudes (DAs) [[l]-[3|] on the previously analysed 24 3 and 16 3 ensembles, 
for which the results can be considered finalised. These will be presented in a forthcoming RBC/ 
UKQCD paper [§J. Preliminary results for the finer 32 3 ensemble will then be presented, with 
a continuum extrapolation in mind. Finally a first look at the use of partially twisted boundary 
conditions for the extraction of distribution amplitudes will be presented. 

The calculations were performed on lattices of sizes 16 3 x 32 x 16 and 24 3 x 64 x 16 with a 
common lattice spacing a" 1 = 1.733(25)GeV and 32 3 x 64 x 16 with cT 1 = 2.2856(28)GeV. The 
gauge field ensembles are drawn from joint RBC/UKQCD datasets using Nf = 2 + 1 flavours of 
domain wall fermions (DWF) and an IWASAKI gauge action. The last dimension, of length 16, 
for each lattice is the extra, 5 , dimension required to implement DWF. The light quark masses 
range from 0.005 to 0.03 (pion masses 331 MeV to 672 MeV) with a fixed strange quark mass, 
am s = 0.04, for the 24 3 ensemble. For the 32 3 ensemble, the light quark mass ranges from 0.004 
to 0.006 (pion masses 290 MeV to 420 MeV) with a unitary am s = 0.03 and partially quenched 
am s = 0.025 [|, [|. 



Ensemble 


am q 


Range 


^Vrneas 


?src locations 


Smearing 


16 3 x 32 


0.005 


900-4480 


180 


0,32,16 


HL-HL 




0.01 


800-3940 


315 


0,32 


GL-GL 




0.02 


1800-3580 


90 


0,32 


HL-HL 




0.03 


1260-3040 


90 


0,32 


HL-HL 


24 3 x 32 


0.01 


500-3990 


350 


0,8,16,24 


GL-GL 




0.02 


500-3990 


350 


0,8,16,24 


GL-GL 




0.03 


4030-7600 


358 


0,16 


GL-GL 


32 3 x 32 


0.004 


760-2410 


86 


0,16,32,48 


LL-LL 




0.006 


500-3240 


82 


0,16,32,48 


LL-LL 




0.008 


500-2960 


61 


0,16,32,48 


LL-LL 



Table 1: Parameters for the 3 datasets. Af meas is the number of measurements for each source position 
f src . The total number of measurements is then N meas X N SIC . In the smearing column XY-XY denotes 
the contraction of two quark propagators with X-type smearing at source and Y-type smearing at sink: 
G=Gaussian wavefunction, H=Hydrogen wavefunction, L=point. 



2. Meson Distribution Amplitudes 

Distribution amplitudes are introduced in the QCD description of hard exclusive processes. 
They encode the non-perturbative QCD effects that occur from factorisation and are important for 
form factors at large q 2 and also for B -decays. They are universal hadronic properties that do not 
depend on the process itself. 

The leading twist DAs for pseudoscalar and vector mesons are defined via vacuum-to-meson 
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matrix elements of quark-antiquark light-cone operators 

(0\q(z)r P Y5^(z-z)q(-z)\K(p))\ zl=0 = f n (ip p ) jT du e^^^U,/!), (2.1) 

(0\q(z)r„75^(z,-z)q(-z)\p(p;l))\ z2=0 = f p m p p^^ j\u e^P^u^). (2.2) 

The path ordered exponential 3? ensures gauge invariance and the DAs are normalised to one 
when integrated over the momentum fraction u carried by the quarks {u = 1 — u). It is useful to 
parameterise the DAs via their moments, 

(S n )M= f'du^MS^) (2-3) 
Jo 

where ^ = u — u is the difference between the momentum fractions. 

These moments appear in matrix elements of local operators with n derivatives. The bare 
moments (denoted (^ M ) bare where n = 1,2) are accessible on the lattice through ratios of two-point 
correlation functions. Taking the simplest example of the first moment pseudoscalar, for large 
Euclidean times t and T — t 



Large t ip P Pfi j, 

*{p»};v^P) = £ xe .>.*< |A v (f,x)Pt(0)|0) ~pV^ ' ^ 

"{ppvyAhP) - £^^(0^(^x^(0)10) Pa {i 1 { ' 

where A and P are the axial and pseudoscalar quark bilinears. The operator with n derivatives is 
given by 

0{Hi 1 ...H,}(x,t) =q(x,t)Y{n D ^...D ^ } q'(x,t) (2.6) 

where the braces in the subscript indicate symmeterisation of the enclosed indices and subtraction 
of traces. 

It is desirable to choose indices so that operator mixing is kept under control and so that there 
are as few non-zero momentum components as possible [g]. We obtain the first moment from 
R{ p4 j. 4 {t , p) (where the index 4 corresponds to the time direction) with p = 1,2 or 3 and a single 
non-zero component of momentum, \p p \ = 2n/L. Similarly the second moment is extracted from 
/?l^ 4 j. 4 (?,/?) with at least two non-zero components of momentum, taking p, ju = 1,2 or 3 with 
p ^ p and \p p \ = \pfi\ =2k/L. A similar method is used to extract moments for the polarised 
vector mesons. 



3. (I 1 ) Results 

For the first moment of the kaon DA, chiral perturbation theory predicts that the extrapolation 
to the chiral limit is linear in m s — m q Our data clearly shows the expected SU(3) symmetry 
breaking effects for all three ensembles. The bare results are plotted in Fig. 1 as a function of the 
mass of the light quark, where m s a = 0.04 for the 16 3 and 24 3 ensembles and m s a = 0.03 (triangular 
points, Fig. 1) and partially quenched m s a = 0.025 (circular points, Fig. 1) for the 32 3 ensemble. 
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0.005 0.01 0.015 0.02 0.025 0.03 0.035 0.04 




0.005 0.01 0.015 0.02 0.025 0.03 0.035 0.04 



am. — am,. 



Figure 1: Linear extrapolation for the first moment K. The solid vertical lines mark the physical point and 
the dashed lines show the error on this point. The left hand plot is 16 3 and 24 3 and the right hand plot is 32 3 
where the circular points are partially quenched. 



We see a similar linear behaviour for (% 1 }k* and so perform the same extrapolation. For the K* we 
see a hint of a finite volume effect, however, where there are results for both volumes at the same 
light-quark mass, they agree within statistical uncertainties. We extrapolate to the physical points 
a(m s -m q ) = 0.0362(16) (16 3 and 24 3 ) and 0.02816(80) (32 3 ) [§ |] for both (^) K and {^) K *. 
The results are given in Table 2. 

4. (£ 2 ) Results 

For the second moments we have some guidance from chiral perturbation theory; that there 
should be no non-analytic dependence at 1 loop and we should fit linearly in m 2 K We see 
no obvious dependence on the quark mass and so perform a linear extrapolation in m\ at a fixed 
simulated m s . The bare results are plotted in Fig. 2. We evaluate the bare second moment values at 
the chiral limit and present these in Table 2. For the second moments we see no obvious indication 
of finite size effects. 





<^ 2 ) b , are 


(or 


( ^ 2) bare 


(or 8 




<or re 


<ar re 


16 3 
24 3 


0.112(5) 
0.125(7) 


0.0228(14)(11) 
0.02377(71)(100) 


0.112(4) 
0.117(5) 


0.02443(96)(107) 
0.0281(13)(14) 


0.110(6) 
0.118(7) 


0.109(10) 
0.118(7) 


0.107(5) 
0.107(4) 


32 3 


0.17(2) 


0.022(1)(1) 


0.14(1) 


0.034(2)(1) 


0.12(2) 


0.15(3) 


0.12(1) 



Table 2: Results for the bare values of moments of the distribution amplitudes. The errors are statistical and 
(for the first moment) due to the uncertainty in the physical point for the chiral extrapolation. Note that bare 
results at different lattice spacings should not be compared directly (as the renormalisation constants depend 
on the lattice spacing). 



5. Renormalisation of results 



We use the Rome-Southampton RI'/MOM scheme to renormalise the relevant operators non- 
perturbatively [Q, ^, |9|]. Multiplying the bare moments by the renormalisation factors, we obtain 
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0.14 - 

0.13 - 

0.12 - 

8 0.11 - 

%t 0.1 - 

0.09 - 

0.08 - 

0.07 - 

0.06 - 



16 3 x 32 
24 3 x 64 




0.005 0.01 0.015 0.02 0.025 0.03 0.035 0.04 



0.2 
0.18 
0.16 
0.14 
0.12 
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0.08 



32 3 x 64 



0.002 0.004 0.006 0.008 0.01 0.012 

am q + am r( , s 



Figure 2: Linear extrapolation for the second moment for the K meson. The solid vertical lines mark the 
physical poi 
quark mass. 



physical point. The Left hand plot is 16 3 and 24 3 and the right hand plot is 32 3 using the unitary strange 



the renormalised results in MS at /i = 2 GeV. It should be noted that the 16 3 and 24 3 lattices have 
different renormalisation factors to the 32 3 . 

The 16 3 and 24 3 results agree with the 32 3 results within or just beyond la (apart from the 
K* first moment). We estimate a formal discretisation error, of 0(a 2 A 1 Q CD ) corresponding to ~ 4% 
(~ 2%) for the 16 3 and 24 3 (32 3 ) ensembles from the 0(a)— improved DWF action and operators, 
which is included in Table 3. The 32 3 results are preliminary and further work will be done to 
confirm its values and errors. Following this, a continuum extrapolation will be made to better 
estimate any discretisation effects. 







it 1 )* 










« 2 >J 


16 3 
24 3 
32 3 


0.25(1)(2) 
0.28(1)(2) 
0.36(5)(2) 


0.035(2)(2) 
0.036(1 )(2) 
0.034(2)(2) 


0.25(1)(2) 
0.26(1)(2) 
0.30(3)(2) 


0.037(1)(2) 
0.043(2)(3) 
0.052(3)(3) 


0.25(1)(2) 
0.25(2)(1) 
0.26(4)(2) 


0.25(2)(1) 
0.27(1)(2) 
0.32(8)(2) 


0.24(1)(1) 
0.25(2)(1) 
0.26(2)(2) 



Table 3: Results for the renormalised values of the distribution amplitudes in MS at fl = 2 GeV. The first 
error is statistical and the second includes systematic errors from m s , discretisation and renormalisation. 



6. Partially twisted boundary conditions 

We used the Z4PSs4 and Z4PSs3 datasets from the RBC/UKQCD K ->■ % form factor runs, 
where the naming convention is as follows: spin-diluted Z(2) x Z(2) noise source and point sink 
with strange quark mass am s = 0.04 and am s = 0.03 respectively. Both sets have am q = 0.005 
and consist of 1180 measurements JlOt]. The gauge field configurations are generated through 
combining sea quarks obeying periodic boundary conditions with valence quarks that have twisted 



boundary conditions Ql ip . The valence quarks therefore satisfy 

y(x k +L) = e i6k y{x k ), k = 1,2,3, (6.1) 
where y/ is either a strange quark or one of the light quarks. 
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We calculate the correlation functions on these datasets with zero fourier momentum and we 
look at cases where only one of the valence quarks is twisted. Therefore the kaon's momentum 
is induced purely by the twist angle of the valence quark. The twist angle is only along one of 
the spatial directions, which is changed regularly in order to reduce correlations. The twist angles 
used for this analysis are di = 1.600 for the light quark for both strange quark masses and d s = 
2.5087(2.7944) for am s = 0.03(0.04) @. The momentum of the meson is then, p = 0/L. 



0.08 
0.07 
0.06 
0.05 
0.04 
0.03 
0.02 
0.01 




Untwisted dafci 
Partially twisted- data (+0.03) 




0.005 0.01 0.015 0.02 0.025 0.03 0.035 0.04 



am* — am,, 



Figure 3: Comparison of using twisted boundary conditions to induce momentum with just using the Fourier 
momentum alone, for the kaon first moment, 24 3 . The partially twisted data shown on the left hand plot is 
for the 0.04 mass strange quark, for which the data-points have been shifted by (+0.03) in the y-direction. 
The partially twisted data on the right hand plot shows the results for both the 0.04 and 0.03 mass strange 
quark. 

Fig. 3 compares the partially twisted data with the untwisted data presented earlier in this 
paper. We see that the results agree well and that using twisted boundary conditions to extract the 
moments of PDAs (at least for the kaon first moment) is possible. It should be noted that we only 
show the cases where the strange quark is twisted. The signal where we twist the light quark is too 



weak to allow a fit. The ratio from which we extract the first moment, see Eq 2.4 , is proportional 
to the momentum and hence the twist angle. The twist is small for the light quark and leads to a 
poor signal for the first moment. 



7. Summary 

We have computed the lowest two non-vanishing moments of the distribution amplitudes for 
the n, K, K*, p and <p mesons, using non-perturbative renormalisation, for two ensembles with a 
common lattice spacing but with different volumes and also for a third ensemble with a finer lattice 
spacing. The renormalised results are presented in Table 3. We do not see any finite size effects 
within errors. The 32 3 results are preliminary and a further investigation will be done to confirm its 
values and errors. We have shown promising results for the extraction of the kaon's first moment 
using partially twisted boundary conditions. 
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